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Abstract 

We study a moving boundary model of non-conserved interface growth that implements the 
interplay between diffusive matter transport and aggregation kinetics at the interface. Conspicuous 
examples are found in thin film production by chemical vapor deposition and electrochemical 
deposition. The model also incorporates noise terms that account for fluctuations in the diffusive 
and in the attachment processes. A small slope approximation allows us to derive effective interface 
evolution equations (IEE) in which parameters are related to those of the full moving boundary 
problem. In particular, the form of the linear dispersion relation of the IEE changes drastically for 
slow or for instantaneous attachment kinetics. In the former case the IEE takes the form of the 
well-known (noisy) Kuramoto-Sivashinsky equation, showing a morphological instability at short 
times that evolves into kinetic roughening of the Kardar-Parisi-Zhang class. In the instantaneous 
kinetics limit, the IEE combines Mullins-Sekerka linear dispersion relation with a KPZ nonlinearity, 
and we provide a numerical study of the ensuing dynamics. In all cases, the long preasymptotic 
transients can account for the experimental difficulties to observe KPZ scaling. We also compare 
our results with relevant data from experiments and discrete models. 

PACS numbers: 68.35.Ct, 81.10.-h, 64.60.Ht, 81.15.Gh, 81.15.Pq 
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I. INTRODUCTION 

Ever since the beginning of the nineteenth century , diffusion-limited growth has 

attracted the attention of physicists, due to its experimental ubiquity and, partly, because 
it is amenable to continuum descriptions that are sometimes solvable. For instance, elec- 
trochemical deposition (ECD) of metals has been and still is a subject of 
intense study during this time due to its (in principle) experimental simplicity and its many 
technological applications: A deposit grows on the cathode when a potential difference is 
set between two metallic electrodes in a salt solution (generally of Cu, Ag or Zn). Another 
interesting system of a conceptually similar type is chemical vapor deposition (CVD) in 
which a deposit grows from a vapor phase through the incorporation of a reacting species 
which attaches via chemical reactions once it reaches the aggregate. CVD is one of the 
techniques of choice for the fabrication of many microelectronic devices, and is currently 
being the object of intense study [h], [ll|, [l2], partly motivated by its use also in emerging 



fields of Science and Technology, such as Microfluidics 131 1 . The practical relevance of ECD 

n n q 

and CVD [4J, [5|, 19[ perhaps makes them appear as two paradigmatic examples of a larger 
class of growth systems (that will be referred to henceforth as diffusive) in which dynamics 
is a result of the competition between diffusive transport and attachment kinetics at the 
aggregate interface. Given that growth dynamics in these processes is not constrained in 
principle by mass conservation, they provide important examples of non-conserved growth 



Despite the great amount of work devoted to these systems, they still pose important 
challenges to the detailed understanding of the very different structures grown under diverse 
conditions, whose geometries range from fractal to columnar. Partial progress has been 
achieved so far through the study of the time evolution of the aggregate surface and its 
roughness 14], ll5|. In particular, a very successful theoretical framework for such type of 
study has been the use of stochastic growth equations for the interface height. Thus, e.g. 
the celebrated Kardar-Parisi-Zhang (KPZ) equation 

^ = V + vV 2 h + ^(\7hf + V (v,t), (1) 

has been postulated as a universal model of non-conserved rough interface growth. In ([I]), v 
is a positive constant, 7/(r, t) is an uncorrected Gaussian noise representing fluctuations, e.g., 
in a flux of depositing particles, and V is the average surface growth velocity. Many times 



the KPZ equation has been put forward as the description of specific experimental growth 
systems based on symmetry considerations and universality arguments. In view of the fact 



171 ] that very few experiments have been reported which are compatible with the predictions 



of the KPZ equation (two examples in ECD and CVD are provided in [18|, [19J), the main 
drawback of such a theoretical approach 1?| is that, in most cases, this coarse-grained 
level of description does not allow to make a connection between the experimental and the 
theoretical parameters, so that it is difficult to assess the cause of the disagreement between 
theoretical description and experimental observations. Moreover, such an approach does 
not allow to include in a systematic way other potentially important physical mechanisms, 
such as e.g. non-local effects (typical of diffusive systems) or fluctuations related to mass 



transport within the di 
In previous work 20 



ute phase. 



21] we have put forward an argument as to why many experiments 



on non-conserved interface growth, such as by ECD and CVD, rarely reproduce the KPZ 
roughness exponents. Namely, morphological instabilities usually occur in these and many 
other growth systems that induce long crossovers making the asymptotic KPZ behavior 
hard to observe. In this paper, we substantiate further such an approach to the problem of 
non-conserved growth, by constructing a model that incorporates the main constitutive laws 
common to diffusive growth systems, from which an effective stochastic growth equation can 
be explicitly derived from first principles. 

The aim of this work is threefold: First, our basic (stochastic) moving boundary problem 
can be explicitly related to realistic CVD and (simplified) ECD systems. We thus provide 
a novel unified picture of these two growth techniques, that makes explicit their common 
features. Nevertheless, due to the generality of the basic constitutive laws that we assume, 
we expect our model to have implications also for different growth procedures that can 
be described as diffusive in the sense described above. Second, we benefit from the model 
formulation in terms of constitutive laws in order to derive the dependence of coefficients 
of the effective interface equation with physical parameters. This result seems to be new in 
the context of diffusive growth, and will allow us to show that the form (and properties) 
of the effective height equation depends crucially on the efficiency of attachment kinetics. 
Specifically, if the kinetics at the surface is instantaneous, i.e., if the particles aggregate 
with probability close to unity when they arrive at the surface, the system can be described 
by a new equation which is morphologically unstable, but that still provides non-KPZ scale 
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invariance of the interface fluctuations at large scales. This new result reinforces our previous 



conclusions 



20 



2l|, on the experimental irrelevance of KPZ scaling in diffusive growth 
systems. On the other hand, for slow interface kinetics the effective evolution equation is the 
well-known (stochastic) Kuramoto-Sivashinsky equation [221 ] . that displays a qualitatively 



similar dynamics, albeit with a long scale behavior that does fall into the KPZ class 23|, [24 , 



251 ]. Third, we will interpret some results from experiments and discrete models under the 
light of our continuum theory both qualitatively and quantitatively. As (the deterministic 
limit of) our model has been profusely tested in the case of CVD, we will mostly consider 
experiments and models from the ECD context. 

The paper is organized as follows. We describe in Sec. II our moving boundary formu- 
lation of diffusive growth systems, including noise terms related to fluctuations in diffusive 
currents and relaxation events. Section III reports a linear stability analysis of the ensuing 
unified model of ECD and CVD. Using a small slope approximation, we derive in Section 
IV a universal nonlinear stochastic equation for the aggregate surface that is numerically 
studied in the novel case of infinitely fast kinetics. Sec. V is devoted to making a connec- 
tion with several experiments and discrete models on diffusive growth systems. Finally, we 
conclude in Sec. VI with a discussion of our results, which will allow us to suggest a reason- 
ably approximate picture of non-conserved growth. Some technical details are given in the 
appendices. 

II. MOVING BOUNDARY PROBLEM 

As it turns out, our description of diffusive growth systems takes a form whose determin- 
istic limit has been long studied in the context of CVD. Thus, we first review the classic 



constitutive equations of CVD |26|, |27|, |28| , and subsequently consider the effect of noise due 
to the fluctuations related to the different relaxation mechanisms involved. Then, we write 
the equations of ECD growth in a form that unifies this technique with CVD. 

A. Chemical vapor deposition 

A stagnant diffusion layer of infinite vertical extent is assumed to exist above the substrate 
upon which an aggregate will grow, see a sketch in Fig. [TJ This approach implies that the 
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FIG. 1: Schematic representation of a model CVD growth system. Black points represent ag- 
gregating units diffusing in the dilute phase. The different transport mechanisms (bulk diffusion, 
attachment and surface diffusion) are indicated in the figure by their corresponding equations. For 
the definitions of the noise terms (q, p, \) see Sec. IIC. 

length of the stagnant layer (typically of the order of cm) is much larger than the typical 
thickness of the deposit (in the range of microns). The particles within the vapor diffuse 
randomly until they arrive at the surface, react and aggregate to it. The concentration of 
these particles, c(x,z,t) = c(r,t), obeys the diffusion equation 

d t c = DV 2 c. (2) 

In the experiments, the mean concentration at the top of the stagnant layer is chosen to be 
a constant, equal to the initial average concentration c a . 

Besides this, mass is conserved at the aggregate surface, so the local normal velocity at 
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an arbitrary point on the surface is given by 



V n = QDVc ■ n - OVo • J. 



(3) 



where Q is the molar volume of the aggregate and n is the local unit normal, exterior to 
the aggregate. The last equation expresses the fact that growth takes place along the local 
normal direction (usually referred to as conformal growth in the CVD literature) and is due 
to the arrival of particles from the vapor [the first term in Eq. ([3])] and via surface diffusion 
(J s stands for the diffusing particle current over the aggregate surface and V s is the surface 
gradient). 

Moreover, the particle concentration, c, and its gradient at the surface are related through 
the mixed boundary condition 



k D (c - c° eq + Tk 



C(x,t) 



DVc-n 



C(x,t) 



(4) 



where c° eq is the local equilibrium concentration of a flat interface in contact with its vapor, 



and C(x,t) is the local surface height. This equation is closely related to the probability of 
a particle to stick to the surface when it reaches it (see below). 

In summary, Eqs. ([2]) and ([3]) describe diffusive transport in the vapor phase and the way 
that particles attach to the growing aggregate. Let us concentrate on the physical meaning 
of Eq. (BJ. The parameter T is related to temperature [29|, [30] as T = jc® q £l/ (fceT), where 
7 is the surface tension (that will be assumed a constant) and n = (d xx () [1 + (d x Q 2 ] 3 ^ 2 
is the surface curvature. The boundary condition (jlj) can be obtained analytically e.g. 
from kinetic theory by computing the probability distribution for a random walker close 
to a partially absorbing boundary. There, the particles have a sticking probability, s, of 
aggregating irreversibly (i.e., attachment is not deterministic). In such a case 31| 



k 



D 



° pr-1 

2-s 



(5) 



where L m f p is the particle mean free path. Assuming that L m f p is sufficiently small we find 
two limits in Eq. If the sticking probability vanishes (s = 0) then Vc = at the 

boundary, so the aggregate does not grow. On the contrary, if the sticking probability is 
close to unity (provided L m f p is small enough), then kjj takes very large values and equation 
01]) reduces to the well-known Gibbs-Thomson relation 29|, [30(] , which incorporates into the 



equations the fact that concentration is different in regions with different curvature. In 



summary, Eq. (jlj) gives a simple macroscopic interpretation of a microscopic parameter, the 
sticking probability, and allows to quantify the efficiency of the chemical reactions leading 
to species attachment at the interface. 



B. Electrochemical deposition 

In an electrochemical experiment, dynamics is more complex than in the CVD system 
as represented above, due to the existence of two different species subject to transport 



f anions and cations) 



32] and an imposed electric field. For a visual re 
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erence, see 



Fig. 12 
351 ] . qualitatively 



Although more elaborate treatments of these can be performed 
the morphological results are similar to the more simplified description we will be making 
in what follows. The virtues of the latter include an explicit mapping to the CVD system 
and explicit experimental verification. 

Thus, in ECD, mass transport is not only due to diffusion in the dilute phase, but also 
due to electromigration and convection. Let C and A be the concentration of cations and 
anions, respectively, then 

d t C = -V • J c , = V • (D C VC - fi c EC - vC) , (6) 
d t A = -V • J a , = V • (D a \7A + fi a EA - vA) , (7) 

where D CA are, respectively, the cationic and anionic diffusion coefficients, /i Cja are their 
mobilities, and E is the electric field through the cell, which obeys the Poisson equation, 

V ■ E = -V 2 = eN A (z c C - z a A)/e, (8) 

with Na the Avogadro constant, ez c and — ez a being the cationic and anionic charges, re- 
spectively, (p the electric potential, and e the fluid permittivity. The velocity v of the fluid 
obeys the Navier-Stokes equation, although we will assume this velocity to vanish in very 
thin cells [36]. 



Another interesting experimental variable is the electric current density, J, given by 

J = F(z c fi c C + z af i a A)E = a(t)E, (9) 

where a(t) is the apparent electric conductivity. Many experiments exploit the capability 
of tuning several parameters while maintaining J constant (galvanostatic conditions), hence 
the relevance of this parameter. 
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Anode 




Cathode 



FIG. 2: (color online) Schematic representation of a model ECD growth system. Cations migrate 
towards the cathode (lower side) while anions migrate towards the anode (upper side) in an infinite 
cell. The different transport mechanisms (bulk diffusion and drift, cation reduction) are indicated 
in the figure by their corresponding equations. 

In order to better understand the way in which the particles evolve in the cell, we need 
to follow their dynamics when the external electric field is switched on. Thus, the cation 
and anion concentrations are initially constant and uniform across the cell and then, once 
the electric field is applied, anions move towards the anode and cations move towards the 
cathode. Cations reduce at the cathode thus forming an aggregate of neutral particles. On 
the contrary, anions do not aggregate, rather, they merely pile up at the anode, which is 
dissolving at the same rate as the cations aggregate in the cathode. Hence, the number of 
cations remains a constant. 

Mathematically, this mechanism of aggregation can be expressed as a boundary condition 
for the cation concentration. Before introducing such a boundary condition we will simplify 
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the set of diffusion equations ©-© following Refs. [37|, [38(. Let us consider that the deposit 
moves with a given constant velocity V, in such a way that, in the frame of reference co- 
moving with the surface, z = is the position of the mean height and z — > oo represents 
the position of the anode (thus, we are dealing with the case in which the height of the 
aggregate is negligible with respect to the electrode separation). Moreover, we will assume 
that the system is under galvanostatic conditions, namely, that the current density at the 
cathode, J, is maintained constant. Thus, the problem can be separated into two spatial 
regions: far enough from and close to the cathode. 

At distances larger than the typical diffusion length, Id = D/V, the net charge is zero, 
so that z a A = z c C. Hence, multiplying Eq. §6§ by z c fi a , and adding Eq. multiplied by 
z a fi c we have 

d t C = DV 2 C, d t A = DV 2 A. (10) 
In both equations we have used the ambipolar diffusion coefficient, given by 

D = VcDg + \loPc / n x 

Hence, mass transport reduces to a single-variable diffusion equation. It is also important 
that the electroneutrality condition (z a A = z c C) implies that the mean interface velocity is 
equal to the anion migration velocit y, t hat is, V = /U a E' 00 , where is the electric field very 

371 . |38| for further details), and then 



far from the cathode (see Refs. 



J-, — z c J c z a 3 a — , (12) 

t 1 — t c 

where t c = fi c / (fi a + fi c ) , and C a is the initial cation concentration. Finally, we must provide 

an equation to describe cation attachment. As a point of departure we will take a relation 

between the charge transport through the interface and its local properties, given by the 

well-known Butler- Volmer (BV) equation ,l33|,|39| 



J=J Q [ e (l-&VcF/*T _ e -(br l+ r ls )z c F/RT Cc/C ^ ? ^3) 

where Jo is the exchange current density in equilibrium, b is a coefficient which ranges from 
to 1, and estimates the asymmetry of the energy barrier related to the cation reduction 
reaction, and rj = A<f)—A(p eq is the overpotential, from which a surface curvature contribution 
i] s has been singled out, of the form 

z c Fr] s 



RT RT 
10 



K = TVlK, (14) 



where we have defined the parameter T = ^/RT in the ECD context. The first term 
on the right hand side of Eq. ffl3l) is proportional to the rate of the backward reaction, 
X — > X n+ + ne~ , and the second one is proportional to the rate of the forward reaction, 
X n+ + ne~ —>■ X. The factor (the concentration at the surface) is due to the supply of 
cations at the surface. Since the flux of anions through the cathode is zero (because they 



neither react nor aggregate) the electric current density at the aggregate surface is on 
to the cations, and the charge current is proportional to the cation current. Hence 

J = VC ■ n 

1 - tr 



,4 



y due 



40] 



(15) 



C 

This equation, combined with Eq. ([13]) . provides a mixed boundary condition which relates 
the cation concentration at the boundary with its gradient. In order to cast it into a shape 
that recalls the CVD relation, we define 

K D = _|L_e-^/«T (16) 
C% = C>^ T , (17) 



and obtain from ffl3|) . ffT5l) 



D. 



c 



VC- n 



K D {C-CfL). (18) 



l-t c 

The coefficient is related with the sticking probability for cations: if the aggregation 
is very effective (large sticking probability) the overpotential is a large negative quantity 
and then Kp grows exponentially. In addition, the concentration C° decreases. Hence, we 
can approximate Eq. ({TBI) by C ~ C° . In the limit when every particle which arrives at 
the surface sticks irreversibly, the solution cannot supply enough particles, C = 0, and the 
current density takes its maximum value. This value of the current is called limiting current 
density. On the other hand, if the sticking probability is small, the system is always close 
to equilibrium and VC ~ 0, so that the net current is zero. 

Finally, we close the system with an equation for mass conservation at the boundary. 
Note that the local velocity of the aggregate surface is proportional to the flux of particles 
arriving to it, therefore 

V n = -OJ c • n = J, (19) 

where Q is the molar volume, here defined as the ratio of the metal molar mass, M, and the 
aggregate mean density, p. For a flat front, V n — V, hence comparing this equation with 
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Eq. flUD we find 

„ M l-t c 

n = — = (20) 

□ 

lip which has been previously proposed theoretically [40J and experimentally 



a relations 
verified 



thus supporting the hypotheses made in this section. 
Surely, the reader has noticed that the last equations resemble those for CVD. To em- 
phasize this similarity, we define new variables and parameters as 



c — R C C, c a — R c C a} cP eq — R c C® q , hp — , (21) 



where R c = D C /[D(1 — t c )). With these definitions, Eqs. (J2])-(|4j) describe (under the physical 
assumptions made above) the evolution of both diffusive growth systems, CVD and ECD. 

C. The role of fluctuations 

The set of equations presented in the previous section describes the evolution of the mean 
value of the concentration so that, formally, we can track the position of the interface at 
any instant. However, it explicitly ignores the (thermal) fluctuations related to the different 
transport and relaxation mechanisms involved. In order to account for these, we define the 
stochastic functions q, p and x as the fluctuations in the flux of particles in the dilute phase 
(— D'Vc), in the surface-diffusing particle current (J s ), and in the equilibrium concentration 
value at the interface, respectively. We choose these noise terms, q, p and Xi to have zero 
mean value and correlations given by 

(fc(r,*)(&(r , ,0> = Qi^-W-O, (22) 
^(r,t)p,(rV)> = P5 v 5 -^==fi > (23) 

( X (r,t)x(r\t')) = I 1 ==A (24) 
where i, j denote vector components and Q, P, and / will be determined from the equilibrium 



fluctuations following |4l|, |42J . Finally, the factor \/l + (d x () 2 in (|23|) , (1241) ensures that the 
noise strength is independent of the surface orientation. 

Thus, the stochastic moving boundary problem we propose to describe diffusive growth 
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has the form 



dtc = DV 2 c - V ■ q, 
Dd n c = k D {c-c° eq + Tn + x) 



+ qn, 



v n = n 



Dd n c - V s ■ 3 S - q • n - V s ■ p 



lim c(x, z\ t) = c a . 



(25) 
(26) 

(27) 
(28) 



In ( 1271) the surface diffusion term, V s • 3 S , is proportional to the surface diffusion coefficient 
D s and the surface concentration of particles v s \ moreover, this term is related to the local 
surface curvature {29], 

V s ■ J s = BV\k. (29) 

In order to determine the values of the coefficients Q, P and I defined in equations 
(EM]) , we use a local equilibrium hypothesis 



42 



43j. To begin with, let us consider an ideal 



concentration c a of randomly distributed particles. The probability of finding n particles in 
a given volume is given by a Poisson distribution. The mean and variance of this distribution 
are c a , hence the concentration c satisfies 



(30) 



(31) 



((c(r, t) - c a )(c(r', - c a )> = c a 5(r - r') 5(t - t'). 

This equation will allow us to determine Q. First, we write Eq. ([25]) as 

d t {c-c a ) = DV 2 (c-c fl )-V-q. 

Let and qk^ be the Fourier transforms of [c(r, t) — c a ] and q, respectively, 

c kw = J dte-™ J rfre- !k ' r [c(r,t)-c fl ], (32) 
q kw = J dte-^ 1 J dre- ik - r q(r,t). (33) 

Writing Eq. ( 13TI) in momentum-frequency space and comparing with the Fourier transform 
of Eq. fl30l . we find that the spectrum of equilibrium fluctuations is (after integrating out 



(c k c_ k ) 



Q_ 

2D' 



(34) 



Hence, using ( 150]) . we obtain Q = 2Dc a . 
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Similarly, in order to determine J, we just note that the equilibrium distribution of a 
curved interface is given by the Boltzmann distribution 

WC})' 



P({C}) ~ exp 



k B T 



(35) 



Ji being a functional that measures the amount of energy needed to create a perturbation 
C(x,t) about the mean interface height. Moreover, as we are assuming a constant surface 
tension 7, 



y/l + (d x Q 2 - 1 



(36) 



fL/2 

j dx 

J -L/2 

provided the perturbation, (, is small enough. The distribution (1351) leads to a fluctuation 



l / dx(d x () 2 , 

Z J -L/2 



spectrum at equilibrium (for a system with lateral dimension L — > 00) of the form 

k B T 



(C*C-*> 



7/c 2 



(37) 



Introducing the boundary condition (1261) into the equation for the velocity ( 1271) and lineariz- 
ing with respect to (, we obtain the following algebraic equation in Fourier space, 

(,kw — - — . n , r;2 Xkid, [oQ 
iuj + UkdI k A 

where we have neglected the surface diffusion terms. Therefore, integrating (Ck^Ck'^') m k' 
and u' and comparing the ensuing fluctuation spectrum to Eq. (13 7p we find 



2Yk B T _ 2c% 



(39) 



Vtk D ^ k D 

Finally, in order to calculate P we assume that the fluctuations due to each relaxation 
mechanism are independent of one another. In this respect, we take the chemical potential 
difference between the interface and the vapor to be given by 



29 



30j (j, = Sl&i, where 5/5( 



denotes functional derivative. Linearizing the equation for the velocity, and considering only 
the contribution due to surface diffusion, we get 



dtC 



V2 1 
+ VSD 



tt 2 v s D S v72 5H 

, - V a — + VSD, 



(40) 



k B T sr ,J1J k B T s 5( 
where V 2 is the Laplace-Beltrami (surface) Laplacian, and tjsd is a noise term related to 

the fluctuations of the surface diffusion current. To ensure that (135]) is the equilibrium 



431 ] . that here 



distribution, this noise term must satisfy the fluctuation-dissipation theorem 
reads 

(vsd(x, t) VS n(x', 0) = 2n 2 u s D s (-V 2 ) 5(x - x')5(t - f). (41) 

Therefore, comparing Eq. (f40|) to Eq. (p7I) we find that — VtV s ■ p = r] SD and consequently 
P = 2LLz/ Q . 
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III. LINEAR STABILITY ANALYSIS 



Eqs. ( |25|) - (!28|) provide a full description of diffusive growth systems including fluctua- 
tions. They thus generalize the classical model of CVD and can also described (through the 
appropriate mapping, as seen above) simplified ECD systems. However, such a stochastic 
moving boundary problem is very hard to handle for practical purposes. In this section we 
will reformulate it into an integro-differential form that will allow us to derive (perturba- 
tively) an approximate evolution equation for the interface height fluctuation, ((x, t). In 
this respect, we will use a technique based on the Green function theorem which has been 



successfully applied to other similar diffusion problems [4l|, |42[ . For brevity, we show here 
the main results leaving the technical details for the interested reader in the appendices. 
Our point of departure is the integro-differential equation 



oo 



- o-(r,t), 

z'=C 

(42) 

where G is the Green function pertinent to the present diffusive problem. The single Eq. (I42p 
relates the concentration at the boundary with the surface height, and is shown in Appendix 
[A] to be actually equivalent to the full set of equations fp5l) -( T28l) . providing essentially the 
so-called Green representation formula for our system [44j. Unfortunately, Eq. (1421) is still 
highly nonlinear and has also multiplicative noise (through the noise term a, see App. |A|) . 
Notwithstanding, it will allow us to perform a perturbative study in a simpler way. 

First, let us consider solutions of Eq. (j42l that are of the form c = c + Ci, where c stands 
for the part associated with the flat (i.e., r-independent) front solution, and c\ is a small 
perturbation of the same order as the height fluctuation ((x, t). Hence, to lowest order in 
the latter and its derivatives (see App. [B]), 

co = , — -■ (43) 

V + k D 

One remarkable feature of the Green function representation is that, from the knowledge of 
the concentration at the boundary, we can extrapolate the value of the particle concentration 
everywhere. Thus, from Eq. ( 14"2"1) and using (1A7I) we find 

c (z) =c a + (c - c a )e- zy / D . (44) 



This equation has been theoretically obtained and experimentally verified by Leger et al. 
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40] 



We now proceed with the next order of the expansion. At this order we already obtain a 
proper (albeit linear) evolution equation for the interface, which moreover contains all the 
noise terms contributions, that is (see App. O . 

dtCk(t)=WkCk(t)+rik(t), (45) 

where Uk is a function (dispersion relation) of wave-vector k whose form may change with 
the values of the phenomenological parameters (see below), and gives the rate at which a 
periodic perturbation of the flat profile grows (if > 0) or decays (if u k < 0) as a function 
of k. Note that, being linear, Eq. (1451) can be exactly solved. 

It turns out that the behavior of can be most significantly studied as a function of 
the values of the kinetic coefficient kp, as anticipated in Sees. Ill Al and III Bl Specifically, we 
analyze separately the case in which surface kinetics is instantaneous (that is, the sticking 
probability is high) and all other cases in which the attachment rate is finite. 



A. Non-instantaneous surface kinetics (hp < oo) 

For a finite value of the kinetic coefficient kjj, we cannot obtain (even in the zero- noise 
limit) the linear dispersion relation Uk in a closed analytic form, unless we perform a large 
scale (k — > 0) approximation. Thus, we can analyze implicitly the zeros of the function 
defined in (IC5I) . which yield the required form of uj^ as a function of wave- vector [451 ] . In 
the large scale limit we find 46] 

u k = a 2 k 2 — a^k 4 , (46) 



where 



Dk 



a-2 



D 



V 



A, (24 



DknlndoIS. 



V 



Vdp 
D 



(47) 



where A = l — do/lu- The two constants appearing in A are the capillarity length (<i = TQ) 
and the diffusion length (Id = D/V). If A < 0, then k = is the only zero of u>k, and 
since a 2 < all Fourier modes Cfc(^) °f the height fluctuation are stable, since they decay 
exponentially in time within linear approximation. 

On the contrary, of A > 0, then a 2 and a 4 are both positive and there is a band of 

r ( i — \ i 1/2 

unstable modes for all k E (0, k*), with k* = ( 1 — y^-jy 1 ) ■ For these values of the 
wave- vector, (k(t) grows exponentially in time within linear approximation. A maximally 
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unstable mode exists corresponding to the maximum positive value of Uk, whose amplitude 
dominates exponentially all other and leads to the formation of a periodic pattern. Under 
these parameter conditions, the dispersion relation (J35J) is that of the linear Kuramoto- 
Sivashinsky (KS) equation (see Fig. [3D \22 \. 

Although the above linear dispersion relation contains 0(k A ) terms, typical of relaxation 



by surface diffusion 29|, [30|, these originate as higher order contributions in which diffusion 



(D), aggregation (ko), and surface tension (r) become coupled. We can also include proper 
surface diffusion into the analysis, for which we simply have to replace 04 by 04 + BQ(V + 
ko)/V, with B as in ( }29j) . which merely shifts k* closer to zero. In this case, the band of 
unstable modes shrinks, which is consistent with the physical smoothing effect of surface- 
diffusion at short lengt 




FIG. 3: Linear dispersion relations given by Eqs. (|46l) (dashed line) and (148)) (solid line), normalized 
by the growth-rate of the most unstable mode, vs spatial frequency k normalized by k*. Both axes 
are in arbitrary units. 



B. Instantaneous surface kinetics (ko — * 00) 



If the sticking probability is essentially one, as seen above kn — > 00. This fast attachment 
condition occurs in many irreversible growth processes 47J. Following a similar procedure 
(and long wavelength approximation) as the one that led us to the KS dispersion relation 
in the previous section, we now get 



u k = D 



T 2 Q 2 



BQ\ , 3TQVk 2 , t/ 9N 



TQV 
D 



r 2 n 2 Btt 



4 



D 



k 



1/2 
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This expression has several interesting limits. For instance, if we neglect surface tension and 
surface diffusion terms (that is, for T — B — 0), then u>k = V\k\, the well-known dispersion 



relation of the Diffusion Limited Aggregation (DLA) model [47j. In this case, every spatial 
length scale is unstable, the shortest ones (large k values) growing faster than the larger 
ones. Thus, in such a case the aggregate consists of wide branches plenty of small tips. 
Moreover, there is actually no characteristic length scale in the system, hence the aggregate 
has scale invariance (that is, it is self-similar). 

If we only neglect the surface diffusion term and, since TQ = do (the capillarity length) 
is typically in the range — 1CT 6 cm, and D/V = Id (the diffusion length) is close to 
lCT 1 — 10~ 2 cm, we can write 

uo k ~V\k\(l-d l D k 2 ) } (48) 



which is the celebrated Mullins-Sekerka (MS) dispersion relation 48J, |49J (see Fig. [3]), ubiq- 



uitous in growth systems in which diffusive instabilities (induced by shadowing of large 
branches over smaller surface features) compete with relaxation by surface tension. This 



dispersion relation has been experimentally verified in several ECD systems 



50 



51 



and has actually been theoretically proposed before for ECD by Barkey et al. 53] (although 
under non-galvanostatic conditions). 

However, in many diffusive growth systems both surface tension and surface diffusion are 
non- negligible; considering again the physical hypothesis do <C Id and a long- wavelength 
approximation, we get 

Uk = V\k\ [1 - (d l D + BQ/2D)k 2 ] - Bttk 4 . (49) 



Nevertheless, there are e.g. some CVD conditions [19j, |54|, for which the vapor pressure 
in the dilute phase is so low that relaxation by evaporation/condensation is negligible in 
practice. In such a case, the dispersion relation is provided by fj49|) with an effective zero 
value for d . 

Finally, there may be physical situations in which quite analogously to the KS case seen 
above, the last dispersion relations [Eqs. (T4"gj) through f!49p ] show the competition between 
mechanisms which tend to destabilize the interface and other which tend to stabilize it. 
From this competition, a characteristic length-scale arises, A m , which grows exponentially 
faster than the others (A m = 27r//c m , with k m being the value for which uo^ is a positive 
maximum) . 
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In summary, we see that, in reducing the efficiency of attachment from complete (ho — > 
oo) to finite (hp < oo), the symmetry of the dispersion relation changes so that non-local 
terms like odd powers of \k\ are replaced by local (linear) interactions. For instance, —k 2 (k(t) 
is the Fourier transform of the local term d 2 ((x, t), while \k\(k{t) cannot be written as (the 
transform of) any local differential operator acting on ((x, t). This result can be understood 
heuristically: if the sticking probability is small, then particles arriving at the interface do 
not stick to it the first time they reach it, but they can explore other regions of the aggregate. 
This attenuates the non-local shadowing effect mentioned above, so that growth becomes 
only due to the local geometry of the surface. 

Although the kr> — > oo limit is a mathematical idealization, for practical purposes we can 
determine under which conditions it is physically attained. Thus, if we introduce Eq. (J33J) 
describing the concentration field for a flat interface into the boundary condition (126]) . the 
latter takes the form 



o _ 

(50) 



C C eq _ c ~ c 



D/k D D/V ' 

The term D/V is the diffusion length; hence, analogously, we can define D/k£> as a sticking 
length. Physically, this length can be seen as the typical distance traveled by a particle 
between its first arrival at the interface and its final sticking site. We can neglect this length 
scale if the sticking probability is close to unity. On the contrary, if ku — > 0, the distance 
that the particle can explore before attaching is infinite. Therefore, taking Eq. (lo*0"j) into 
account, we can say that the kr> — > oo limit describes accurately the problem when kr> 3> V, 
and in such a case the diffusion length and the capillarity length determine the characteristic 
length-scale of the system. 



IV. NONLINEAR EVOLUTION EQUATION 

In this section we proceed one step further with our perturbative approach by including 
the lowest order nonlinear contributions to Eq. (1451) . If we evaluate the Green function of 
the problem at the boundary we find 



G(r 



r ,t 



e(r) 

A-kDt 



exp 



{x-x'f (C-C + Vrf 



ADt 



ADt 



(51) 



where r = t — t' . Expanding the last term in the argument of the exponential as a series 
in ( we get (£ — ( r ) 2 + 2V(( — C0 r + V 2 t 2 . The second and third terms were already taken 
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into account before in the linear analysis, so the only nonlinear contribution in Eq. (|5"T|) is 
related to the first term, (£ — (') 2 , that introduces a correcting factor exp[— (£ — C') 2 / (4-Dt)] 
which is only significant when £' ~ £, hence we can replace (' — ( by the lowest term in its 
Taylor expansion, to get 

By incorporating this contribution into the formulae of appendices \K\ through (UJ the conse- 
quence can be readily seen to be the addition of a mere term equal to the (Fourier transform 
of) V(d x () 2 /2 to the right hand side of Eq. (j4"51) . resulting into an evolution equation with 
the form 

d t Ck(t) = co k Ck(t) + ^Af[(] k + Vk(t), (53) 

where jV[C]fc is the Fourier transform of jV[C] = (d x () 2 . Note that the nonlinear term 
obtained is precisely the characteristic KPZ nonlinearity as in Eq. (CQ), that appears here with 
a coefficient equal to half the average growth velocity, agreeing with standard mesoscopic 



arguments 



161 ] . Moreover, as is well known 14|, this is the most relevant nonlinear term 



(asymptotically) that can be obtained for a non-conservative growth equation such as Eq. 
( |42l . hence any other nonlinear term will not change the long time, long scale behavior of 
the system. As we have shown above, for small sticking Uk is given by Eq. ( 1461) and for 
large sticking is given by Eq. (j4"5]) . In all cases the noise correlations involve constant terms, 
as well as terms that are proportional to successively higher powers of k. By retaining only 
the lowest order contributions in a long-wavelength and quasistatic approximation, from Eq. 
(1C8|) we obtain 

(■nkuVk'w) = (no + U 2 k 2 )5(k + k')5{u + u% (54) 

where 

' Vc (l + ^), k D <oo 
U ={ V kD ' , (55) 

Vc° eq , k D — ► oo 

n2 = j 2^ 2 c (^ + 4) + + & 2 , k D < oo 

^ 2D 2 c° eq /V + 2D s v s , k D ^oo 

where Eq. ( 1431) for Co is to be used in the case of finite kinetics. In general, parameter 
n provides the strength of non-conserved noise, while n 2 measures the contribution of 
conserved noise 55| to the interface fluctuations. The presence of conserved noise can 
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and 



naturally modify some short length and time scales of the system but, in the presence of 
non-conserved noise, it is known to be irrelevant to the large scale behavior. Thus, H2 will 
be neglected in the numerical study performed below. 

For the case of finite kinetic coefficient, the evolution equation (|53p is the stochastic 



generalization of the KS equation [23|, |24], |25j. In the case of fast attachment kjj — > 00, 
the resulting interface equation combines the linear dispersion relation of MS with the KPZ 
nonlinearity. In this sense it employs two "ingredients" that seem ubiquitous in growth 
systems, so we find it remarkable the fact that (to the best of our knowledge) its detailed 
dynamics has not been reported so far. The purpose of the next section is to report a 
numerical study of this equation in order to clarify the similarities and differences to its 
finite attachment counterpart . 



A. Numerical results 



1. The pseudo- spectral method 



As noted above, although the shape of the nonlinear Eq. (|53|) is common to both stick- 
ing limits, their dispersion relations make them very different physically. Thus, while the 
linear terms of the KS equation are local in space, linear terms corresponding to MS disper- 
sion relation cannot be written in terms of local spatial derivatives. Therefore, we cannot 
perform a standard finite difference discretization in order to implement its numerical in- 
tegration 55j. Rather, in order to integrate numerically Eq. (1531) we will resort to the 



so-called pseudo-spectral methods that make use both of real and Fourier space representa- 
tions. Such techniques have been successfully used e.g. in many instances of the Physics of 



Fluids 



561 ] , a nd are bein g u sed more recently in the study of stochastic partial differential 



equations [57 



58 



59|, 



60 



61]- 



As we are interested in the qualitative scaling properties of Eq. (|53|) . rather than, say, in 
a quantitative comparison to a specific physical system, we introduce positive constants z/, 
K, B, and A, that allow us to write the equations in the more general form for each limit: 



dtCk(t) = (uk 2 -Kk 4 )C k + ^Af[Q k + Vk (t), 
d t ( k (t) = (v\k\-K\k\ 3 )( k + ^Af[C} k + Vk(t), 



(57) 
(58) 
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where, in order to stress the similarities between the two interface equations, we have ne- 
glected in (|58p the 0(k A ) terms, given that their stabilizing role is already played by the 
0(|A;| 3 ) term. In what follows, we will only refer to the new parameters, which will later be 
estimated in the next section when we analyze different experimental conditions. 

In order to integrate efficiently Eq. (I5"5j) we use a pseudo-spectral scheme. This numerical 



method is detailed in 



6l| and employs an auxiliary change of variable that allows to estimate 



the updated value of (k as 

Ck(t + At) = e^ At Uk(t) + -AtAr[C(t)], 
where the noise term r k {t) is conveniently expressed as 



+ r k {t), (59) 



e 2w fe At _ 1 

r k (t) = x \ (n + W) v k (t), (60) 

with Vk{t) being the Fourier transform of a set of Gaussian random numbers with zero mean 
and unit variance Regarding the explicit calculation of the nonlinear term in Eq. (I59I) . 
we perform the inverse Fourier transform of —ik( k and take the square of it in real space, 
so that we explicitly avoid nonlinear discretization issues 57j. However, in this procedure 



aliasing issues arise [56( , that we avoid by extending the number of Fourier modes involved 
in the integration, and using zero-padding, see details in 0, Q ■ 



2. Numerical integration 

As mentioned, the noisy KS equation, Eq. (|571) . has been extensively studied in the 
literature, so we will concentrate in this section on the novel Eq. ( 158]) . Its behavior will help 
us to understand the evolution of the interface in the ko —>■ oo cases. 

The linear regime is very similar for both equations, as one might naively expect. In 
fact, we find that both systems feature similar power spectral densities (or surface structure 
factors), S(k,t) = (Cfc(^)C— , see Fig. HI This can be easily understood by inspection of 
the analytic result for S(k,t) obtained from the exact solution of Eq. ( ]45l) . 

p 2uj k t _ i 

S(k,t) = U — , (61) 

where the corresponding dispersion relations are displayed in Fig. [3j By simple inspection 
of Fig. (jlj) one is tempted to say that both Eqs. fl57|) and fl58l have a similar behavior so 
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FIG. 4: (color online) Power spectral density, S(k,t), from numerical simulations for a L = 1024 
system with v = K = 1, Ho = 1CP 2 , A = (i.e., linearized equations), averaged over 10 3 realiza- 
tions, at times t = 4, 8, 10 and 20 for: (a) Eq. (f5T|) (circles) and (b) Eq. (|58|) (circles). Blue solid 
lines represent the exact solution, Eq. (|6ip . for each case. All axes are in arbitrary units. 

that, a priori, KPZ scaling might be expected in the asymptotic regime also for the latter. 
However, for this fast kinetics equation, we show in Fig. [5] that the growth exponent j3 
characterizing the power-law growth of the surface roughness or global width W(t) ~ t 13 , 



where W 2 (t) = (1/L 2 ) £ fc S(M) [U, [15|, is much larger at long times than the expected 
KPZ value @kpz = 1/3. A rationale for such a long-time behavior of Eq. (158!) can be 
already provided by simple dimensional analysis. Thus, under the scale transformation 
t b z t, k -> b- l k, Cfc -> b a+1 ( k , Eq. becomes 

dtCk = b z - l v\k\( k - V-*K\k\X k + b a+z - 2 ^M[Q k + V' 2 -*- 1 ' 2 ^. (62) 

If we introduce the exact one dimensional KPZ exponents (a = 1/2, z = 3/2), it is easy to see 
that in the hydrodynamic limit (that is, when b — > oo) the most relevant term in the equation 
is not the KPZ term but, rather, the lowest order linear term \k\( k . Preliminary dynamic 
renormalization group calculations 62j seem to provide the same result. The present scaling 
argument provides moreover the exponent values a = /3 = z = l&t the stationary state. 
These values are compatible with those obtained from numerical simulations of Eq. (1581) . 
as displayed in Figs. and The numerical values we obtain for the exponents are a = 
1.00 ± 0.05, (3 = 1.05 ± 0.05 and z = 0.95 ± 0.05, thus they are in good agreement with the 
ones predicted by dimensional analysis. In order to check the consistency of our numerical 



estimates 
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641 ]. in the inset of Fig. [6] we show the collapse of the power spectrum density 
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using these exponent values. Collapses are satisfactory, including the behavior of the scaling 
function for the width, indicated by a solid line in the inset of Fig. [51 The discrepancies 
in the collapsed curves for large kt l l z values is due to the existence of a short scale scaling 
different from the asymptotic one. 

From our numerical results, we conclude that the non-linear regime as described by Eq. 
(158|) . corresponding to instantaneous kinetics, is very different from that for slow kinetics, 
as represented by the noisy KS equation. Thus, in both cases the KPZ nonlinearity is able 
to stabilize the system and induce power law growth of the surface roughness, associated 
with kinetic roughening properties. However, the universality class of Eq. ( 1581) is not that 
of the KPZ equation but, rather, it is a new class completely determined by the |fc| term in 
the linear dispersion relation. Note moreover that the new exponents associated with the 
asymptotic regime for this equation fulfill accidentally the Galilean scaling relation a + z = 2 



14 



16( 1 . We describe this property as accidental due to the fact that Eq. (1581) is not Galilean 
invariant. The easiest way to confirm this is to check for the scaling behavior of a stable 
version of Eq. (1581) in which we take negative v values for which all modes are linearly 
stable. Fig. [7] shows on the left panel the power spectral density in this case, that behaves as 
S(k) ~ 1/k for long distances. Moreover, on the right panel of this figure we plot the time 
evolution of the global roughness for the same system, that behaves as W(t) ~ logt for long 
times before saturation. From these data we conclude the exponent values are z = 1, and 
a = (log), which do not satisfy the scaling relation implied by Galilean invariance. Such 
a stabilized version of Eq. ( 158]) has been studied in the context of diffusion-limited erosion 
65|. 



V. COMPARISON WITH EXPERIMENTS AND DISCRETE MODELS 

In this section, we focus on the applications of the model equations to understand and 
explain some experimental results, qualitatively and quantitatively and, besides, compare 
our results from continuum theory to relevant discrete models of diffusion limited growth. 
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FIG. 5: (color online) Global width vs time for a system with = i^T = A = 1, IIo = 10 -2 obtained 
numerically for Eq. (|58|) . for increasing system sizes, L = 32, 64, 128, 256, 512 and 1024, bottom 
to top, averaged over 10 3 realizations. The red dashed line is a guide to the eye with slope 1.05 
suggesting the asymptotic value of j3. Inset: Collapse of W(t) using a = 1.00, j3 = 1.05, and 
z = 0.95. The blue line has slope 1.00, showing the consistency of our estimate for a. Axes in the 
main panel and in the i 



FIG. 6: (color online) Power spectrum vs spatial frequency k with the same parameters as in Fig. 
[5] and L = 1024. Different curves stand for different times (the one at the bottom is for the earliest 
time). The dashed line is a guide to the eye with slope —3, compatible with a = 1. Inset: Collapse 
of S(k,t) using a = 1.00, (3 = 1.05, and z = 0.95. Axes in the main panel and in the inset are all 
in arbitrary units. 
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FIG. 7: (color online) Left: Power spectral density vs wave vector for Eq. (|58|) with v = —1, K = 1, 
A = 4, Ho = 10, L = 8192, and averages over 100 noise realizations. The dashed line is a guide 
to the eye with slope —1 [that is, a = (log)]. Right: Global roughness vs time for the same 
system as in the left panel. Note the semilog representation. The dashed line is a guide to the eye 
representing W(t) ~ logi before saturation. All axes are in arbitrary units. 

A. Experiments 

From the experimental point of view, it is very difficult to determine if the dispersion 
relation characterizing a specific physical system is the MS or, rather, the KS one, because 
they are both very similar except very near k = (see Fig. [3]). In principle, such a dis- 
tinction would be very informative, since it could provide a method to assess whether the 
dynamics is diffusion limited (kr> — » oo), or else reaction limited (ko <C V). Several previous 



theoretical works in this field 
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67J predict, qualitatively, a KS dispersion relation, while 
other 53(, predict a dispersion relation similar to the one of MS. All those studies are not 
necessarily incompatible with one another because they are considering different experimen- 
tal conditions. Moreover, to our knowledge, there are only a few experimental reports in 
which the dispersion relation is measured, and in most cases the authors presume that it can 



be accounted for by MS 



51 



521 ] . Notwithstanding, within experimental uncertainties 
that are specially severe at small wave-vectors, they could all equally have been described 
by the KS dispersion relation. 

Another way to distinguish which is the correct effective interface equation that describes 
a given system could be the value of the characteristic length associated with the most 
unstable mode that can be measured (which in fact would be essentially the characteristic 
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length-scale that could observed macroscopically for the system). Unfortunately, from Eq. 
f|4"T|) one has that, for d <C Id, then k m = (2Zu(io)~ 1 ^ 2 for the KS and k m = (SZdcZo) -1 / 2 for 
the MS dispersion relation. Hence, except for a constant numerical value of order unity, both 
cases provide a characteristic length scale that depends equally on physical parameters, while 
the significant parameter, hp, only modifies u m , namely, the characteristic time at which the 
instability appears (which is about uj^ 1 ). Thus, in order to clarify the nature of the growth 
regime (diffusion or reaction limited), one should rather study the long time behavior of the 
interface. 



Despite these difficulties, we can sil 



in the literature. Leger et al. 



38 



40 



try to interpret some experimental results reported 



68| have presented several exhaustive works dealing 



with ECD of Cu under galvanostatic conditions. In addition to properties related to the 
aggregate, they also provide detailed information about the cation concentration. Their 
main result in this respect is that such concentration obeys experimentally Eq. (PHI) as we 
anticipated above. It can also be seen that the product aggregates are branched and that 
the topmost site at every position x defines a rough front growing with a constant velocity. 
In Fig. 8 of Ref. 40j the authors plot the branch width against the diffusion length. From 
that figure, it is clear that A m and Id are not linearly related. Moreover, they seem to better 

1/2 

agree with the present prediction from either of our effective interface equations, A m oc , 



than with the linear behavior argued for in 



40]. 



Another important experimental feature is the relation between A m and C a . From Eq. 
201) we know that Q = (1 — t c )/C a , then do = FIT = Q'-f/RT, so that the characteristic 
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length scale, A m , will be proportional to C a . Thus, one would expect the branches to be 
narrower as we increase the initial concentration, consistent with the patterns obtained by 
Leger et al. 

We will try also to interpret the ECD experiments of Schilardi et al. [jij. They have 
measured the interface global width considering that the topmost heights of the branches 
provide a well defined front. Thus, the time evolution of the global width, or roughness, 
presents three different well defined regimes: A short initial transient, which cannot be 
accurately characterized by any power-law due to the lack of measured points, is followed 
by an unstable transient. We consider the system unstable in the sense that the average 
interface velocity is not a constant but, rather, grows with time. Finally, the system reaches a 
regime characterized by exponent values that are compatible with those of KPZ the equation, 
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while the aggregate grows at a constant velocity. These regimes again resemble qualitatively 
the behavior expected for the noisy KS equation fl57|) . 

Moreover, we can check whether the order of magnitude of the experimental parameters 
is compatible with our predictions. First, we estimate k m from mean width of the branches 
within the unstable regime. This width is about 0.05 mm 18|, hence k m ~ 1.3 x 10 3 
cm -1 . Besides this, the mean aggregate velocity at long times is V ~ 2 x 10 -4 cms -1 , 
and the diffusion coefficient is D ~ 10 -5 cm 2 s" 1 , so that the diffusion length is about 
Id = D/V ~ 0.05 cm. These magnitudes allow us to calculate the capillarity length d = 
l/2l£)k% n ~ 5 x 10~ 6 cm (which is of the same order as we considered in our approximations 
above, and much smaller than the diffusion length). Furthermore, the instability appears at 
times of order l/uj m . In the experiment this time is about 6 min. Thus, u m ~ 3 x 10 -3 s _1 . 
Finally, as u m = then ko ~ 6 x 10~ 8 cms -1 <C V, which provides a consistency 

criterion for the validity of our approximations and of our predictions. 

As a final example, let us perform some comparison with the mentioned experiments of 
"q^ . At short times, they obtain compact aggregates with exponent 



Pastor and Rubio 
values 
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52| a = 1.3 ± 0.2, ai oc = 0.9 ± 0.1, z = 3.2 ± 0.3, and /3 = 0.4 ± 0.08. This 
means that the interface is superrough (a > 1). After this superrough regime, the aggregate 
becomes unstable and the dispersion relation has the MS form. The fact that the aggregates 
are compact (and not ramified) seems to show that surface diffusion is an important growth 
mechanism in these experiments (which is also consistent with the small value of the velocity, 
V ~ 4 /mi/mm). Consequently, we can use Eq. ( 158]) with an additional surface diffusion 
contribution [— Bk 4 (k{t)] to understand the behavior of the experiments. We choose the 
parameters v = 0.25, K = U = 1 and A = n 2 = (because the velocity is small and we are 
only interested in the short time regime), for several values of B between and 1 in order 
to determine the influence of surface diffusion in the growth exponents. Other parameters 
only change the characteristic length and time scales of the experiment. The exponents 
thus obtained numerically are (with B = 0.75) (3 ~ 0.39 ± 0.02, and a ~ 1.3 ± 0.1, which 
are (within error bars) equal to the experimental ones. As we have pointed out above, this 
superrough regime is followed by an unstable transient characterized by the MS dispersion 
relation, as has been also observed in other ECD experiments by de Bruyn 5(|, and Kahanda 



et al. 



5l|- 



More recently, additional ECD experiments have been reported under galvanostatic con- 
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ditions. E.g. in Ref. {(| three growth regimes can be distinguished: a first one at short 
times, in which a Mullins-Sekerka like instability is reported, is followed by a regime in 
which anomalous scaling (namely, the roughness exponents measured from the global and 



local surface widths differ, a jk a ]n 



15 



63 



64j takes place, and finally at long times or- 



dinary Family- Vicsek scaling 14| is recovered. Similar transitions to and from anomalous 
scaling behavior have been also reported in Ref. jg]. Our present theory does not predict 
the anomalous scaling regimes reported by these authors. This may be due to the small 
slope condition employed in the derivation of equations (158]) and (!57|) . However, we want 
to emphasize two points in this respect: (i) In Ref. js] the authors report a transition from 
rough interface behavior to mound formation. These mounds can be obtained by numerical 
integration of equation ([SSD in 2+1 dimensions [to] . (ii) As we will show in the next section, 
the theory is in good agreement with a discrete model of growth in which anomalous scal- 
ing is clearly reproduced. Hopefully, a numerical integration of the full moving boundary 
problem (j42~l) would capture the anomalous scaling regime, and this will be the subject of 
further work. 



B. Discrete Models 



As mentioned above, kr> is related with the sticking probability through Eq. (jSJ). This 



probability acts as a noise reduction parameter Tlj] in discrete growth models, as was shown 
e.g. in 0, S] for the Multiparticle Biased Diffusion Limited Aggregation (MBDLA) 
model, used to study ECD growth. In particular, MBDLA has been seen to describe quanti- 



tatively the morphologies obtained in [181 ] . In MBDLA, by reducing the sticking probability 
the asymptotic KPZ scaling is indeed more readily achieved, reducing the importance of 
pre-asymptotic unstable transients, as illustrated by Fig. 8 of [741 ] . Hence, noise reduction is 
not a mere computational tool for discrete models but, rather, it can be intimately connected 
with the surface kinetics via Eqs. (j"5J) and (T261) . 

MBDLA with surface diffusion also predicts the existence of a characteristic branch width. 
This is shown in Fig. [8] in which the power spectral density is plotted and compared with 
the one obtained from the noisy KS equation, proving the equivalence between both descrip- 
tions of ECD. These results are reinforced by the fact that, as shown in Ref. [74J], the cation 
concentration obeys Eq. (14*41 . and the branch width dependence on the cation concentration 
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FIG. 8: (color online) Power spectrum obtained from MBDLA simulations (solid line) with r = 0.5, 



741 ]. by permission), and Eq. (|57p with u = 1, K = 1/4, 



p = 0.5, s = 1, c = 1 (taken from Fig. 15 in 
A = 40, Ilo = 10 -2 , L = 512 (dashed line), averaged over 10 5 realizations. For the sake of clarity, 
the latter has been vertically offset. Straight blue lines are guides to the eye having slope —4. Axes 
are in arbitrary units. 

— 1/2 

is consistent with the relation A m oc C a 1 . Moreover, in simulations of MBDLA, an un- 
stable transient was found before the KPZ scaling regime, being characterized by intrinsic 
anomalous scaling, as recently observed in the experimental works by Huo and Schwarzacher 
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76] . As mentioned, probably the absence of such an anomalous scaling transient in our 



continuum model is related to the small slope approximation, and we expect to retrieve it 
from a numerical integration of the full moving boundary problem. 



VI. CONCLUSIONS 



In this paper we have provided a derivation of stochastic interface equations from the 
basic constitutive laws that apply to growth system in which aggregating units are subject 
to diffusive transport, and attach only after (possibly) finite reaction kinetics. Our derivation 
seems to be new for this class of systems, and allows to relate the coefficients in the effective 
interface equation with physical parameters, like the sticking parameter, physical surface 
tension and size of aggregating units, etc. We have seen that the shape of the equation 
describing the time evolution of the aggregate interface changes as a function of the sticking 
probability. For very high interface kinetics non-local shadowing effects occur, while for 
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finite kinetics non-local shadowing yields to morphological instabilities of a local nature. 
Thus, qualitatively the behavior of the system for generic parameter conditions roughly 
consists of an initial transient associated with morphological instabilities in which typical 
length scales are selected (thus breaking scale invariance), that is followed by a late time 
regime in which the interface displays kinetic roughening. However, the universality class of 
the latter differs, being of the KPZ class only for finite attachment kinetics, while it becomes 
of a new, different, non-KPZ type for infinitely fast attachment, as predicted by the new 
interface equation we obtain in the latter condition. While in previous reports [l5, 1?], 0] 
we have interpreted the long unstable transients as a potential cause for the experimental 
difficulty in observing KPZ scaling, our present results go one step further in the sense that 
for fast attachment conditions we do not even expect KPZ universality in the asymptotic 
state, due to the irrelevance of the KPZ nonlinearity as compared with the \k\(k{t) term in 
Eq. (EHD. 

Comparison of our continuum model with experiments and discrete models seem to sup- 
port the above conclusions. Nevertheless, our results are in principle constrained by a 
small-slope approximation. In view (specially under the fast kinetics conditions) of the 
large roughness exponent values that characterize the long time interfaces as described by 
our effective interface equations, it is natural to question whether the same scenario holds 
for the full (stochastic) moving boundary problem. Moreover, our small slope equations 
do not account for anomalous scaling, which is otherwise seen in experiments and in the 
MBDLA model, so that integration of the complete system ( T25l) -( l28l) seems indeed in or- 
der. Technically, systems of this type pose severe difficulties even to numerical simulations 
(mostly related to front tracking in the face of overhang formation). Thus, one needs to 
rephrase the original continuum description fl25l) - fl28l) into an equivalent formulation that is 
more amenable to efficient numerical simulation, such as e.g. a phase-field model 77|. We 
are currently pursuing such type of approach, and expect to report on it soon. 
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APPENDIX A: THE GREEN FUNCTION TECHNIQUE 



This technique has been used in other similar problems, such as solidification 41] or 
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epitaxial growth on vicinal surfaces 4^|, and is based on the use of the Green theorem 
to transform an integral extended over a certain domain (in our case, say, the region between 
the electrodes (see Fig. [9]) to an integral that is evaluated precisely at the moving boundary 
(the aggregate surface). 

Let us consider that the distance separating the electrodes and the lateral size of the 
system are both infinite, so that the only part of the dashed line in Fig. [9]whose contribution 
is non vanishing is the aggregate surface. The Green function related to Eq. (1251) is the 



A 




FIG. 9: Integration domain, A, and its boundary (solid line) used in Green theorem. The infinites- 
imal arc length along the moving boundary ((x,t) is given by ds. 

solution of 

(Jp + DV 2 - V^j G(r -r',t- t>) = -8(v - r')5(t - t'), (Al) 

where we have made a change of coordinates to a frame of reference moving with the average 
growth velocity, V. To evaluate G, we use its Fourier transform 

G ( r _ r ',t _ t ') = J dwe^-V j d'ke^-^G^, (A2) 
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where r = ix + zi and k = k x x + k z z. Hence, Eq. (lAlj) becomes 

Glr, 



Dk 2 + icu -iVk z ' 
with k 2 = k 2 + k 2 . Therefore, integrating Eq. (1A2[) we find 



(A3) 



G(r - r', * - t') = %■ exp 



{x-x'f (z-z' + V{t-t'))' 



AD(t - V) 



W(t-t>) 



(A4) 



Q(t — t') being the Heaviside step function. In the following, we will use for brevity r = t — t'. 
It can be straightforwardly seen that the following relations are satisfied: 



lim G(r — r'; r) = S(r — r') 

T— »0+ 



lim G(r — r'; r) 



0. 



T— > — OO 



dt' / dx'G 

OO J — OO 



exp [-(2 - z>)V/D\/V si 2 > z', 
1/V si z < z', 



(A5) 
(A6) 

(A7) 



We can now rewrite Eq. ff25l) in terms of the variables r' and t', 

£-nV- V °L) </ = -v>.4 (A8) 

Adding Eq. flA8} multiplied by G(r -r',t- tf) to Eq. (jAl]) multiplied by d = c(x' , z' ; t') , 
and integrating t' in (—00, t — e] and r' in the set A = (—00, 00) x [((x', t'), 00), see Fig. [91 



t-e 



8 /"*" e f 8 



t-e 



.4 



dt' / dV^c'G) -V I dt 1 I dV^c'G) + 



df I ( /V(rV-6'- G'V'V 



t-e 

dt 7 / d 2 r'GV ■ q 



2//-ir7/ „/ 



(A9) 



The first term on the left hand side can be easily evaluated with the use of Eqs. ( 1A5I) and 
( fA6l) . Thus, 



lim 

e^O 



t-e 



/ d 2 r'—(c'G) = c(r;t) + / dt' / dx' ^ 



dt' L 



c'G 



(A10) 



where £' stand for £(x',t'). Analogously, we can integrate the second term on the left hand 
side of Eq. (IS9jl . using fjATjl 



t-e 



iim-y / dt' / dV^-fc'G) = -c + y / dt' / dx' 



c'G 



z'=(> 



(All) 
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Finally, using the identity dV 2 G — GV 2 d = V • (c'V'G — GV'c'), and applying the Green 
theorem on the domain A, we get 



MmD f 6 dt' [ d 2 r' (c'V 2 G - GV 2 d) = -D f dt' [ d 

J-oo J A ' J -co JC 



s> (c>^ - G— 
dn' dn' 



(A12) 



ds' being the arc length, see Fig. [91 from which we obtain the integro-differential equation 



c(r,t) 



dt' 



dx'[V+ ^pl c'G 



D, d.Ac> d S-G d ^ + ^ 



dn 1 



dn' 



z'=C 



where 



(A13) 



(A14) 



dt' / dx I dz'GV' ■ q' 
> J-oo Jc' 

is a term related to the diffusion noise. 

As we seek to determine c everywhere, we need to know its value at the boundary. 
Considering the limit in which r belongs to that boundary (hereafter, we will denote it as 
rj), the term dG/dn' in Eq. (1 A 13ft is singular. Fortunately, this singularity is integrable, as 
a result of which we obtain an additional term as l42|l 



,dG 



+ 



,dG 



(A15) 



This leads to Eq. (f42|) of the main text, where we have omitted the subindex b for notational 
simplicity. 

APPENDIX B: ZEROTH ORDER CALCULATION 



Writing c = c + C\ in Eq. (|42j) we get 



Co ci _ 
2 + 2 ~ ° a ~ 



CO J —CO 



dt' I dx 



z'=C 

(Bl) 



Note that, at this order, ds' = y 1 + (d x Q 2 dx' ~ dx' . We also linearly expand G, so that 



G(r-r',t-t r ) ~ (l-V 



where 



AuDt 



cxp 



2D 

(x - x') 2 V 2 r 
4Dt 4D 



(B2) 



(B3) 
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In order to determine the concentration at the boundary, we must write our main equation 
in terms of cq and c\. Then, with the use of the boundary condition ff26l) we find 

<9(c + ci) 



D- 



dn 



k D (c + cx- c° + Td 2 x C + X) + q • n. 



(B4) 



Finally, putting Eqs. (|B2]) and flB3} into flE} we find 



2 + 2 ~ Ca 



L Co ^-(A; D +y/2)( c [ ) + c ' 1 )+A ;i? y^(c-C / )+^c° ? (i- 



2D 



K - 



r , i C-Cfv 2 l 



c 



a(r,t), (B5) 



z'=C 



with a new noise term 



/t /»oo r /*oo 

dt' / cfo' (k DX ' + q' • n')G° + / d2'G°V • q' 
-OO J— OO L ■/ f 



(B6) 



Despite the apparent complexity of these new equations, Eq. flBlj) is linear so that Fourier 
transforming it we get the following algebraic equation which relates all the zeroth order 
terms 



C a + 



k D c, 



eq 



V 



8{k)8{u) - 



V 



(P ku being the Fourier transformation of Co and 



ADui + AD 2 k 2 + V 



-1/2 



(B7) 



(B8) 



which can be easily inverted yielding Eq. fj43l . 



APPENDIX C: FIRST ORDER CALCULATION 



From the results obtained in Apps. [A] and [B] we can find an evolution equation which 
relates c\ w and Cfcw Thus, 



■'kuj 



- tuc - (k D + + -L (k D c - k D c° eq ) (-^- -v)+ Tk D k 2 + 



-V*)- -Dk 



1 

akuj being the Fourier transformation of a(r,t), hence 
1 V 



C ku 



2G L 



+ k 



D 



)- 






4) 




QD 


K 2G L 





+ k D Tk 2 



Ckw — °ku 



(C2) 
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This equation has two unknowns; hence, in order to solve it, we also need to expand Eq. 

(127)) in powers of C\ and (. Thus, 

1 ( ^ 2 B A ik 

c ku, = 7^j— + Tfc + — k (ku, - Xku + -j—Pkuj- (C3) 
H2fc D k D J k D 



Combining both Eqs. jC2j and ()C3j) . we find 
with 

The new noise term is the projection of all the noise terms onto the boundary, and is 
given by 

(I V \ ( ik \ 

= DAi + J Xk . ~ ^pUD^J + k D ) + £ dz> (ikfa + g z k Al ] ) exp(-At ) z'),(C6) 



where 



a(±) - 1 + V (C7) 
u ~ 2DG° kw 2D' iUj 

These equations allow us to calculate the noise correlations (in Fourier space) 

</WW> = [^lA^I 2 + 2Dc ( ^ + l fPA + 2D s u^(d^\A^ + fc 2 , + 2DJfc D Re( 
fcD \ 2Re(A^J) J k D \ 

x5{k + k')5{uj + uj'). (C8) 

Note that, in principle, Eq. flC4j) provides us with all the information of the system at linear 
order, namely, the power spectrum of £ (by evaluating T k ~J and then integrating out the 
temporal frequency u) or the height-height correlations (integrating out the spatial frequency 
k). 

In order to gain insight about the implications of this expansion in (, we write Eq. ( 1C4j) 
as the Fourier transformation of a Langevin equation for the interface height, 

[iuj - u^Cfco, = Tjku, (C9) 

uj k being a function of k which we must specify from T kLU , and rj kLJ being a noise term related 
with Pku, that is also delta-correlated, 

fatoW) = n(fc)5(fc + k')S{u + u'), (CIO) 
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where II depends, in general, on k and provides the magnitude of the noise for each Fourier 
mode, see Sect. IIVI Finally, by Fourier transforming the time frequency, we find the linear 
stochastic partial differential equation (H5]l . 
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